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BACKWARD ASSEMBLY PLANNING WITH DFA ANALYSIS 
AWARDS ABSTRACT 

An assembly planning system that operates based on a recursive 
decomposition of assembly into subassemblies, and analyzes 
assembly cost in terms of stability, directionality, and 
manipulability to guide the generation of preferred assembly 
plans. The planning in this system incorporates the special 
processes, such as cleaning, testing, labeling, etc. that must 
occur during the assembly, and handles nonreversible as well as 
reversible assembly tasks through backward assembly planning. In 
order to increase the planning efficiency, the system avoids the 
analysis of decompositions that do not correspond to feasible 
assembly tasks. This is achieved by grouping and merging those 
parts that can not be decomposable at the current stage of 
backward assembly planning due to the requirement of special 
processes and the constraint of interconnection feasibility. The 
invention includes methods of evaluating assembly cost in terms 
of the number of fixtures (or holding devices) and reorientations 
required for assembly, through the analysis of stability, 
directionality, and manipulability. All these factors are used 
in defining cost and heuristic functions for an AO* search for 
an optimal plan. 
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5 ANALYSIS 

BACKGROUND OF THE INVENTION 
Origin of the Invention: 

The invention described herein was made in the performance of work 
under a NASA contract, and is subject to the provisions of Public Law 
96-517 (35 USC 202) in which the contractor has elected not to retain 
title. 

Technical Field: 

15 The invention relates to a method and apparatus for performing back- 

ward assembly planning and operation using backward search of assem- 
bly sequences. 

Background Art: 

20 Assembly planning can be based either on the forward search of as- 

sembly sequences or on the backward search of assembly sequences, given 
geometric and topological information on the product to be assembled. 
The backward search of assembly sequences is often implemented in the 
form of disassembly planning, assuming that an assembly sequence can 
25 be obtained by reversing the order of the corresponding disassembly se- 
quence. Although there are cases where assembly sequences may not 
be the reverse of the corresponding disassembly sequences, disassembly 
planning is widely used in the research community due to its advantage 
in planning efficiency: The disassemblability of a part or a subassembly 
3 Q directly implies the satisfaction of precedence relationship, whereas, in 
the forward search, the satisfaction of precedence relationship between a 
pair of mating parts may not be known immediately until an exhaustive 
search is completed. 

Much effort has been expended to make assembly planning or disas- 
sembly planning more autonomous, more efficient, and closer to reality. 
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Early work on interactive assembly planning was concerned about for- 
mulating a necessary and sufficient set of questions to be answered by 
the human designer that leads to the complete specification of precedence 
relationships with the minimum number of question- and- answering oper- 
ations. Subsequently, the emphasis has been moved to automatic reason- 
ing of geometric interference or path planning as a means of identifying 
precedence relationships in assembly. The feasibility of automating ge- 
ometric reasoning in assembly planning is fueled by the development of 
powerful computer-based geometric modeling and reasoning tools. The 
automatic test for the existence of a mating path in disassembly or as- 
sembly is still computationally very expensive; and this significantly lim- 
its the number of parts that can be handled by an automatic assembly 
planner, since the number of feasible assembly sequences may grow expo- 
nentially, so does the number of required tests for path existence. One of 
15 the current research goals is thus to achieve planning efficiency by avoid- 
ing unnecessary tests for path existence, as proposed by such methods 
as the logical inference of path existence based on the previous results 
on path planning, the test of necessary conditions for successful mating, 
which are simple to implement, prior to the application of path planing, 
20 the use of path finding heuristics. 

It becomes clear that, to bring assembly planning closer to reality, we 
need to consider not only identifying the geometric interference in part 
mating, but also many other factors which affect assembly planning in 
a real factory environment: special processes such as testing, labeling, 
25 painting, grinding, etc. to be performed during assembly; physical con- 
straints due to magnetic and electrostatic interference; stability of parts 
and subassemblies which directly affect assembly cost due to required 
fixtures; part and subassembly manipulability; the required number of 
reorientations during assembly; tolerance propagation; as well as the 
30 physical constraints imposed by assembly layouts. The above factors 
should eventually be incorporated into the process of identifying prece- 
dence relationships as well as selecting preferred assembly sequences. 
For instance, the required special processes, the physical interference, 
and the stability requirement can be used in precedence indentification 
together with the geometric interference condition. While, the stability, 
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the manipulability, the assembly direction and reorientation, and the tol- 
erance associated with parts, subassemblies, and their matings can be 
used in selecting prefered assembly sequences. However, the selection of 
a cost-effective assembly sequence is not simple and remains as one of 
major research issues in assembly planning. This is not only due to the 
fact that evaluating the above mentioned conditions is difficult and com- 
putationally expensive but also the optimality in assembly planning may 
not be achievable through the local optimization of the above conditions 
but requires the consideration of global effect. 

The present invention includes a process for backward assembly plan- 
ning and performance which takes the above issues into consideration. 


15 SUMMARY OF THE INVENTION 

A method of selecting and performing a sequence of sub- assembly 
operations to construct an assembly begins by identifying all direct sub- 
assemblies of the assembly by identifying potential subassemblies of the 
assembly, and, for each potential subassembly, determining whether the 
20 potential subassembly is a direct subassembly of the assembly by: (1) de- 
termining whether there is a common axis of separation between the po- 
tential subassembly and the remainder of the assembly, ( 2 ) determining 
whether the potential subassembly is stable, and (3) determining whether 
there is a direct unobstructed path between the potential subassembly 
25 and the remainder of the assembly. The identifying potential subassem- 
blies is performed by identifying all liaisons connecting pairs of nodes of 
each potential subassembly of the assembly, and, for each of the liaisons, 
merging liaisons for which there is no accessibility for interconnection 
between the subassembly and the assembly and merging liaisons which 
30 have no force-deliverable paths to the assembly. It is further performed 
by identifying special processes of the assembly of the type comprising 
testing and cleaning which require certain nodes of the assembly to be 
processed together during a special process. The identifying potential 
subassemblies further consists of grouping together nodes of the assem- 
bly which must be processed together during special processes of the 
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type comprising cleaning and testing of subassemblies, whereby to pre- 
vent identification of potential subassemblies comprising corresponding 
nodes of the grouping. The overall process further consists of identifying 
precedent constraints corresponding to required sequences of assembly 
^ of parts of the assembly, and the identifying potential subassemblies fur- 
ther consists of grouping together nodes of the assembly which must 
be assembled in a hierarchical sequence corresponding to the precedent 
constraints, whereby to prevent identification of potential subassemblies 
comprising corresponding nodes of the hierarchical sequence. The deter- 
^ mining whether potential subassembly is stable consists of determining 
that one of (a) the potential subassembly has more only one transla- 
tional degree of freedom and only one rotational degree of freedom, the 
one rotational degree of freedom being about the translational axis only, 
or (b) any floating of disconnected liaison of the potential subassem- 
15 bly includes a node of the remainder of the assembly. The determining 
whether the potential subassembly is stable consists of determining that 
the potential subassembly is unstable if either the potential subassembly 
is unstable or the remainder of the assembly is unstable. The process 
further consists of computing a cost of mating and assemblying each 
20 direct subassembly, by determining the number of part reorientations 
required to construct the subassembly and determining the number of 
mating operations required to construct the subassembly. The process 
further consists of representing the assembly as a tree of AND and OR 
nodes beginning with all direct subassemblies as bottom nodes of the tree 
25 and concluding with the assembly as the top node of the tree, the tree 
comprising branches connecting adjacent ones of the nodes, the branches 
corresponding to alternative decompositions of the assembly and com- 
puting a local cost of each OR node of the tree and eliminating branches 
corresponding to OR nodes of higher costs, and performing an assembly 
30 operation of subassemblies corresponding to the tree after the eliminat- 
ing. The local cost computation consists of determining the relative cost 
due to the stabilization of the direct subassembly and the remainder of 
the assembly by using holding devices, and determining the relative cost 
due to reorientation operations. 

BRIEF DESCRIPTION OF THE DRAWINGS 
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Fig. 1 is an exploded view of a flashlight assembly employed as a tu- 
torial example of an application of the present invention. 

Fig. 2 is a liaison graph of the flashlight assembly of Fig. 1. 

Fig. 3 illustrates the special process forest of the flashlight assembly 
of Fig. 1. 

Fig. 4 illustrates two assemblies in which an assembly sequence is not 
the reverse of a disassembly sequence. 

Fig. 5 illustrates the clustering process of the flashlight assembly of 

Fig. 1. 

15 Fig. 6 illustrates the special process forest of Fig. 3 after removing 

the root node. 

Fig. 7 illustrates how the internal freedom of motion of a subassembly 
is determined. 

20 

Fig. 8 illustrates an AND-OR tree employed in the invention as the 
search space for the AO* algorithm. 

Fig. 9 illustates the AND-OR tree corresponding to the flashlight as- 

25 sembly of Fig. 1. 

Fig. 10 illustrates the merging and grouping operations corresponding 
to the liaison graph of Fig. 2. 

30 Fig. 11 illustrates the overall process of the invention. 

Fig. 12 illustrates a system for carrying out the invention. 


DETAILED DESCRIPTION OF THE INVENTION 
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1 Assembly Representation 

An assembly A is a cluster of parts assembled together by a certain 
assembly sequence, which maintains a particular geometric relationship 
5 among parts. More specifically, a cluster of parts, P, is considered as an 
assembly, iff 

1) Every part in P has at least one geometric constraint in relation with 
other parts in P. 

2) There exists at least one assembly sequence which can put all the 
parts of P together to satisfy the geometric constraints imposed on 
themselves. 

3) P is stable in a sense that it can maintain its geometric constraints 
either by itself or through the aid of holding devices. 

The above definition extends the notion of conventional assembly to 
include a disconnected and non-self-stable assembly that may be gener- 
ated during the assembly of the product. 

Formally, an assembly A is represented by 4-t.uples: 

A = {P(A), Gl{A), Gp(A), n(A)}, where P{A) represents the set of 
parts constituting A, Gi(A ) the attributed liaison graph representation 
of A, Gp(A ) the special process and constraint forest associated with A, 
and n(A) the set of all the feasible assembly sequences for A. 

1.1 Attributed Liaison Graph 

A liaison is said to exist between a pair of parts if one part constrains 
the freedom of motion of the other by a direct contact. A liaison graph 
is a graph, G, G = ( N,E ), with N representing a set of nodes, and E 
representing a set of edges. A node n, n £ N , is assigned to each part 
of an assembly, and an edge e, e E P, is assigned to a liaison between a 
pair of parts. An attributed liaison graph is a liaison graph with frames 
attached to individual nodes and edges of the liaison graph to describe 
the attribute associated with a node or an edge. A part frame attached 
to each node describes the attributes associated with a part, including 
1) the part geometry and coordinate, 2) the mating volumes and the 
contact surfaces as part features, and 3) the physical properties of a part 
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30 
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such as weight. An edge frame attached to each liaison describes at- 
tributes associated with the liaison. The attributes of a liaison consist 
of 1) the parts in associated with the liaison, 2) the mating features and 
mating directions, 3) the mating type such as insertion and place-on, 
and 4) the interconnection type such as Attachment, Force-fit, Connec- 
tors/Retainers, Push-and-Twist, Screw, Glue, or Welding. Thus, G L (A) 
contains information on the topology of part configurations, the geom- 
etry and relative pose of parts in A, the interconnection mechanisms of 
part connections, and the local freedom of motion in part mating. 


1.2 Special Process and Constraint Forest 


The assembly of a product involves not only the interconnection of parts 
to form required liaisons but also the execution of special processes such 
as testing, cleaning, surface treatment, painting and packaging during as- 
sembly. In addition, there may exist assembly constraints that should be 
satisfied in order to prevent electrostatic, electromagnetic, and thermal 
interference, as well as mechanical vibrations and chemical pollutions, 
during assembly. The accomplishment of special processes and the satis- 
faction of assembly constraints impose precedence relationships in assem- 
bly. Especially, assembly precedence relationships should be generated 
in such a way as to guarantee the execution of multiple special piocesses 
in a certain precedence order. 

Special processes and assembly constraints associated with -4 are rep- 
resented by the special process and constraint forest, Gj,(-4), consisting 
of the special process forest, G S (A), and the assembly constraint forest, 


G C (A). 

The special process forest, G S (A), 
following properties: 


consists of a set of trees having the 


1) A node n- of G S (A) represents a special process, 5/. n- is associated 

with a tuple (i?,£f), where Pf represents a set of parts 
involved in 5,-, and £• represents the union of the parts involved 
in Si and all the special processes that should precede Sr. £• = 
iYU{Uj£j},Vnj: n ) is a c b.ild of n- (a recursive definition). 

2 ) A branch 6- connecting the two nodes, n* and nj, represents the 


precedence relationship between Si and Sy. Sj < Si if nj is a child 
of nj. Special processes corresponding to the sibling nodes of a tree 
or to the nodes of different trees have no precedence relationship. 

The assembly constraint forest, G C (A), consists of a set of trees having 
the following properties: 

1) A node nj of G C (A) represents a liaison Z, of Gl{A). n f is associated 

with a tuple where Pf represents a pair of parts involved 

in Z,, E* represents the union of the parts involved in Z, and all the 
liaisons that should precede /,•: E* = i^U{Uj E^}< Vnj, n c j is a child 
of nj. 

2) A branch bj connecting nj and n c j represents a precedence relationship 

between Z, and lj". lj < Z,- if nj is a child of nj. Liaisons corresponding 
to the sibling nodes of a tree or to the nodes of different trees have 
no precedence relationship. 

Example: Flashlight Assembly 

The flashlight assembly shown in Fig. 1 consists of a body(BD), 
two dry cells(Bl,B2),a bulb carrier(BC), a bulb(BL), a receptacle(R), 
lens(G), and a cover(C). Fig. 2 shows the attributed liaison graph repre- 
sentation of the flashlight assembly. Fig. 3 illustrates the special process 
forest associated with the flashlight assembly. It shows that the cleaning 
of the receptacle(R) should precede the testing of the bulb unit(BC,BL), 
such that BC, BL, and R should be clustered in a subassembly during 
backward assembly planning. 

2 Backward Assembly Planning 
2.1 Definition 

A subassembly of an assembly A , S,| A, is an assembly that can be gen- 
erated in one of the feasible assembly sequences for A. A direct sub- 
assembly of A, Sf\A, is defined as a subassembly of A, which can be 
directly assembled with A - Sf\A at the last step of one of the feasible 
assembly sequences for A. The assembly planning of an assembly A, 
A = {P(A),G l (A), Gp(A), n(A)}, is the process of generating a set 
of assembly sequences fla(A) C 11(A) based on the given P(A), Gi{A) 


and Gp(A ), and the criteria for selecting preferred assembly sequences. 
The backward assembly planning of A , BAP(A), represents a particular 
method for achieving the assembly planning of A , based on the recursive 
identification and selection of preferred direct subassemblies. BAP(A) 
first identifies and selects a direct subassembly of A , Sf |A(or a set of 
direct subassemblies of A, {5f|A,z = 1, • ■ • , m}) and decomposes A into 
Sf\A and A - Sf\A (or {Sf\A and A - S?\A,i = Then, 

it recursively applies the process of decomposition to the subassemblies 
generated by the previous decomposition until no further decomposition 
can be applied. 

2.2 Backward Assembly Planning vs. Disassembly Planning 

BAP(A) differs from disassembly planning in that BAP (A) can handle 
the case where an assembly sequence can not be obtained from the re- 
verse of a disassembly sequence. For instance, a sequence of operations 
to disconnect a liaison of force-fit interconnection type is often quite dif- 
ferent from the reverse of a sequence of operations to interconnect the 
liaison. As shown in Fig. 4a, the snaps of part A should be widened in 
order to disassemble part B from part A, which may require a new tool. 
A more interesting example is given in Fig. 4b, where we can have the 
following disassembly sequence: C-A-D-B, since Screw C can be disas- 
sembled before Cover A is removed. However, the reverse of the above 
disassembly sequence, B-D-A- C, can not be an assembly sequence, since 
the inability of holding part D during assembly of screw C prohibits such 
an assembly sequence. BAP (A) is able to identify that the direct sub- 
assembly of (A,B,C,D) is A but not C, and generate a feasible backward 
assembly sequence, A-C-D-B. 

2.3 Identification of S’flA 

The core of BAP(A) is the identification of Sf|A. For the identification 
of Sf\A, we introduce the following conditions that a cluster of parts 
P\A, P\A C A, can be a direct subassembly of A: 

1) Accessibility Condition : P\A is accessible by a gripper or a tool. 


2) Stability Condition : P\A and A - P\A are stable, i.e., they either 
are self-stable or can be stabilized by external devices. 

3) Local Mating Motion Condition : All liaisons between -P|-4 and .4 — 
P\A have at least one common axis of separation. 

4) Path Existence Condition : P\A can be brought to .4 — -P|-4 from 
the free space for mating. 

5) Interconnection Feasibility Condition : P\A can be interconnected 
to A — P\A by applying the interconnection operations defined for 
the liaisons between P\A and A — P\A. 

6) Process Constraint Condition : jP|A meets the constraints defined 
by the special process and constraint forest of A, Gp(A). 

To avoid the unnecessary evaluation of computationally expensive con- 
ditions, such as the path existence condition, we hierarchically evaluate 
the above conditions. The following sections describe the details of the 
method for testing the above conditions. 

3 Part Clustering 

The identification of Sf\A can be accomplished by checking the above 
conditions for individual P|A’s defined by the cut-sets of Gi(A). Since 
the number of cut-sets of G i(A) is often large and grows exponentially 
with respect to the number of nodes in Gl{A), testing all the cut -sets of 
Gl(A) results in inefficiency in assembly planning (especially due to the 
high cost of testing the path existence condition). 

The number of the cut-sets that require the test of path existence 
condition may be reduced considerably if we first select those cut-sets 
that satisfy the interconnection feasibility condition, the process con- 
straint condition, and the stability condition, prior to the evaluation of 
the path existence condition. For the selection of those cut-sets satisfy 
the interconnection feasibility and process constraint conditions, we first 
construct an abstract liaison graph, Gl(A), from Gl{A), as follows: 

1) Merging those parts of A that cannot be interconnected, although 
it is assumed that those parts are already brought into their mating 
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position. This process identifies the liaisons of Gi{A) that violate 

the interconnection feasibility condition, by reasoning whether the 

preconditions for a particular interconnection mechanism associated 

with a liaison can be satisfied. 

5 

2) Grouping those parts that should belong to the same subassembly 
to meet the requirement of special processes. 

The parts merged together become a supernode, whereas the parts 
grouped together become a group node in Gl(A). Then, the cut-sets 
of Gl{A), where a super node and a group node are considered as a 
single node, satisfy the interconnection feasibility and process constraint 
conditions. 

3.1 Interconnection Feasibility and Part Merging 

15 A liaison between P\A and A — P\A cannot be completed, even if P\A 
can be brought into its mating position, if the following conditions are 
not satisfied: 

1) The feasibility of applying tools and installing connectors required for 

the interconnection of L. 

20 

2) The feasibility of applying a force, while maintaining stability, re- 

quired for the interconnection of 

The test of the first condition can be done by checking the existence 
of an open channel to the designated part locations, through which tools 
25 and connectors can be operated without geometric interference. 

The test of the second condition requires reasoning on the force deliv- 
ery to a liaison /, through intermediate liaisons. To be more specific, let 
us first categorize a liaison into one of the following three classes: 

3Q 1) A liaison is said to be floating if there exists no physical force hold- 
ing the parts (associated with the liaison) together. For instance, a 
liaison with the interconnection type of “attachment” is a floating 
liaison. A floating liaison may or may not be stable, depending on 
the geometric constraints imposed on the parts associated with the 
liaison. 
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2 ) A liaison is said to be rigid, if there exists physical force holding the 
associated parts together, by which the liaison becomes self- stable 
even under the presence of external force. For instance, a liaison with 
the interconnection type of “force-fit”, “welding”, or “connectors” 
may be classified as a rigid liaison. 

3) A liaison is said to be firm if there exists a physical force holding the 
associated parts together, by which the liaison becomes self-stable, 
although it may exhibit a deformation or a freedom of motion under 

IQ the presence of external force. 

A liaison is associated with its local freedom of motion: 

The local freedom of motion of a liaison U, LFM(lj). connecting the 
two parts, Pi and P 2 , is represented by the freedom of motion of P\ 
against P 2 , XPM(/,; P 1 IP 2 ) or the freedom of motion of P 2 against Pi, 
15 LFM{U\ P 2 |Pi), where LFM(l { ; Pi|P 2 ) or LFM(U\ P 2 |P] ) is represented 
in terms of the coordinates of the assembly to which Pi and P 2 belong. 

A local freedom of motion of a liaison may be changed after the in- 
terconnection is completed: A rigid liaison completely loses its local 
freedom of motion after the interconnection is completed. A firm liaison 
20 may show a local freedom of motion, when the force applied to the liaison 
is more than a certain threshold. To distinguish the local freedom of mo- 
tion of U after theinterconnection from that before the interconnection, 
LFM + (k ) is used to represent LFM{U) after the interconnection. Note 
that, due to the orthogonality between motion space and force space, 
25 a liaison connot deliver force to the direction where a local freedom of 
motion exists. 

To test the feasibility of applying force required for the interconnection 
of , let us define the following: 

Definition: An Accessible Node , A-node 
30 

A part is accessible by a tool if it is reachable and graspable by a 
tool for an assembly operation. A node is accessible if any of the 
parts forming the node is accessible. 

Definition: An Access Path to a node n, A-path 
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An access path to a node n, A-path, is a path starting from an A- 
node and ending with the node n without having any other A-node 
in the middle of the path. By definition, an A-node has an A-path 
to itself. An access path is represented by an ordered set of liaisons 
or parts on the path. A node may have one or more A-paths. A 
pair of A-paths, A-path(ni) and A-path(ri2), n\ ^ 71-2, are said to be 
independent each other, if they share no common part. 

Definition: The Internal Motion Space, Al[A-path(n)], and the Static 
Force Space , J r [A-path(n)], associated with A-path(?i) 

The internal motion space, Ad[A-path(n)], of A- pat.h( 7 i) is defined 
by the union of the local freedom of motion of individual liaisons in A- 
path(n), and represents the flexibility that the configuration of parts 
along A-path(n) can be deformed by an external force, with the first 
part (corresponding to A- node) and the last part (corresponding to 
the node n) fixed in space. Note that At [ A-path ( 11 )] is a function of 
the external force applied to A-path(n) when A-path (71 ) includes firm 
liaisons, since the external force determines which firm liaisons can be 
broken. As will be explained shortly, the external force to be applied 
to A-path(n) is given as the force required for the interconnection of 
the liaison which n is associated with, and is subject to the evaluation 
of force-deliverability. 

Assume that A-path(ri) is represented by an ordered set of liaisons, 
{h,h, ' ’ ' Jr}, with /,• formed by a pair of nodes, (77,,, 77 ;,), and that 
(71^,77,2) is ordered along A-path(n) in the direction toward n. Then, 
A 4 [A-path(n)] =U- = i LFM(W, n.-Jn^,). The static force space, F[A- 
path(n)], of A- path(Ti) defines the static force that can be deliv- 
ered to the node n through A-path(n). A[A-path( 7 i)] is represented 
by the orthogonal complement of Af[A- path( 7 ?)]. i.e., A[A-path( 77 )] 
={±x,±y, ±z,±ip, ±9,± 0} — \j r i=1 LFM(li; n^\n t .,). 

Definition: A Force-Deliverable A-path 

A-path(njj) is said to be force-deliverable to n i} for the liaison U of 
(77^,71,2), if A[A-path(?7 11 )] includes the force required for the inter- 
connection of 77*2 tO 77 ,, . 
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The interconnection feasibility condition of a liaison of (n,-, , w* 2 ) is now 
transformed to the verification of the existence of an independent force- 
deliverable A-path for n,j and n,- 2 . The force-deliverability of an A-path 
to rii l or rq 2 depends on the force required for the interconnection of /,-. 
That is, M [A-path(nq )] or M [A- path(n, 2 )] is determined by the freedom 
of motion of individual liaisons along the A-path under the presence of 
an external force equivalent to the force required for the interconnection 
of U. For instance, in case U is floating, the amount of force required by U 
is negligible, and the decision of the internal motion space of the A-path 
^ is based solely on the floating liaisons along the path. Note that, in this 
case, the force-deliverability of /, becomes equivalent to the feasibility of 
maintaining stability during interconnection. 

As a summary, we present the following merging principle: 

15 Merging Principle : 

The liaison /, of (n tl , n,- 2 ) can not be listed as a cut-set, and conse- 
quently Ui x and rii 2 should be merged together, if one of the following 
conditions is true: 

1) It is not feasible for the tools and connectors (required for the 

20 interconnection of /,) to access the designated locations. 

2) Either n,, or n, 2 has no independent force-deliverable A-path, 
including the case where either n,-, or n, 2 has no A-path at all. 

The following algorithm summarizes the process for part merging: 

25 Algorithm: Part Merging 

Input: Gl(A). 

Output: Gi{A) with super nodes. 

Method: Step 1 . Put all the liaisons of Gl(A) in Open set. Iden- 
3 q tify A-nodes of Gi{A). 

Step 2. If Open set is empty, then stop. 

Step 3. Select a liaison U from Open set in an order starting 
from an A-node, and remove Li from Open set. 

Step 4. Check whether l t requires tools or connectors to com- 
plete the interconnection. If not, go to Step 6. 
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Step 5. Check the accessibility of tools or connectors to the des- 
ignated locations. If not, merge n,i and n ,2 associated with /, 
and go to Step 2. 

Step 6. Check whether nn and have independent, force-deliverab 
A-paths. If not, merge n,i and n,- 2 and go to Step 2. 

Example : Part Merging of the Flashlight Assembly 
The A-nodes of the flashlight assembly are identified as BD and C 
(refer to Fig. 2). To show how each liaison in Fig. 2 can be tested for 
10 the merging, let us consider lg as an example: 

For lg of (C,R), A-path for C is {C}, whereas A-paths for R are {C,R}, 
{C,G,R}, {BD,B2,BC,R}, {BD,B1,B2,BC,R} and {BD,B1,B2,BC,BL,R}. 
lg does not require a tool or a connector to complete the interconnection. 
The part, C, has a force-deliverable A-path, {C}; however the part, R, 
does not have a force-deliverable A-path independent of {C}. This is 
because all A-paths for.R from BD can not deliver force to R due to the 
floating liaisons luhihihi and h • Therefore, lg is merged. 

Fig. 5 shows the result of the merging process. 

3.2 Special Process Constraints and Part Grouping 

20 

The recursive determination of Sf\A in backward assembly planning 
should support the execution of special processes in an order specified 
by G S (A) and satisfy the assembly constraints specified by This 

requires the identification of those parts of A that should be grouped 
together, or should not be divided into Sf|A and A — 5f|A, in order to 
^ implement the special process forest G S (A) and the assembly constraint 
forest G C {A) as follows: 

1) Given G S (A), we first determine which special processes should be 
executed with A, i.e., the selection of one or more special processes 
2 Q executable with A from G S (A) in a top-down order from the root 

nodes (note that there can be multiple trees in G S (A)). We then re- 
move out those nodes corresponding to the selected special processes 
from G S {A), resulting in G S (A). Since the special processes remain- 
ing in G' S {A) need to be accomplished in the later stage of backward 
assembly planning, we should group those parts in the accumulated 
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part list of each of the root nodes of G' S (A). 

2) In order to preserve the precedence relationship defined by G c { A), we 
should group those parts in the accumulated part list of each child 
^ of individual root nodes of G C (A). Note that a liaison represented 

by a root node of G C (A) is eligible for decomposition at the current 
stage of backward assembly planning, and is exempt from grouping. 

Grouping Principle : 

For a given A, Gl(A), G s (A ) and G C (A ), we group those nodes of 
10 Gl(A) that belong to either of the following lists: 

1) Eo associated with each root node, 71 q., of G' S (A), where G S (A) 
is obtained by removing out from G S (A ) those nodes that are 
selected for processing with A. 

2 ) Ei, associated with each child, n\., of individual root nodes of 

15 Gc(A). 

Example: Part Grouping of the Flashlight Assembly 
First, G' S (A) is obtained in Fig. 6 , since the testing of the whole 
assembly should be done with A. 

20 Then, {B1,B2} and {BC,BL,R}, the accumulated part lists of indi- 

vidual root nodes of G' S (A), become groups and {BC,BL.R} becomes 
another group, as shown in Fig. 5. {G} is a single part and is not indeed 
in the group list. 
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4 Identification of Direct Subassemblies 

The problem of finding direct subassemblies is now transformed into the 
problem of finding valid cut-sets of the abstract liaison graph. G z. ( A ) . 
A cut-set, 7 ,-, decomposes Gl{A) into disjoint subgraphs, G ^ [A.~,j) and 
G£ 2 (A, 7f), where Gj (A, 7 ,) or G£ 2 (A, 7«) may not be a connected graph, 
but may be a collection of multiple connected subgraphs, G^^A^',) = 
{G^ 1; (A, ji), j = 1,2, ■••,/!} or G£ 2 (A, 7 i )={Gf 2 ; (A, 7 1 ),i = l,2,---,/ 2 }- 
A cut-set, 7 i, is valid if: 1) Each connected subgraph. G/ J (A, 7 ,:), j = 
1,2, •••,/!, G?'(A,7i), j = 1,2- generated by the cut-set, 7 , in- 
cludes at least one A-node. 2) There exists a collision free path between 
the cluster of parts, P\A, corresponding to ££^4,7,) and 
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A — P\A, corresponding to Gi(A,ji) for their mating. 

The first condition ensures that P\A and A — P\A, can be handled by 
a tool for assembly, whereas the second condition ensures that Sf\A can 
be placed in its mating position. For the test of the second condition, 
let us define the following: 

1) The predicate LM(P\A, A — P\A ) is true if all the liaisons between 
P\A and A — P\ A has at least one common axis of separation. 

2) The predicate PE(P\A, A — P\A) is true if there exists a path along 
which the cluster of parts P\A can be brought to its mating position 
without colliding with the rest of the assembly. 

Since the premise that LM(P\A, A — P\A ) is not true is a sufficient 
condition that PE(P\A, A — P\A ) is not true, the test of LM{P\A, A — 
P\A ) preceding the test of PE(P\A, A — P\A) can provide a considerable 
reduction in the number of the costly PE(P\A, A — P\A) test. 

Algorithm: Identification of Direct Subassemblies 

Input: Gl(A), A-Node Set = {all the A-nodes of GT(A)}, Tested Cut- 
set List— {0}, Valid Cut-set List={0}. 

Output: A list of direct subassemblies specified in Valid Cut-set List. 

Method: 

Step 1. If A-node Set is empty, stop. 

Step 2. Select an A-node from A-node Set, and remove it from A- 
node Set. 

Step 3. If there exists a cut-set 7,-, 7* £ Tested Cut-set List, such 
that G S i (A, 7 i) includes the A-node selected in Step 2 and G L -{A, 7,) 
includes at least one A- node of G L (A ), then continue. Otherwise, 
go to Step 1. 

Step 4 . Put 7 i into Tested Cut-set List. 

Step 5. Test 7 f for LM. If it is false, go to Step 3. 

Step 6 . Test 7,- for stability. That is, if one of the subassemblies 
produced by 7, is unstable or the subassemblies produced by 7, 
have no common assembly pose among themselves, go to Step 3. 
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Step 7. Test 7 ,• for PE. If it is true, put 7, in Valid Cut-set List. 
Otherwise, go to Step 3. 

The details of the method for testing the stability of 7, will be pre- 
5 sented in the next section. 

5 Subassembly Evaluation 

The direct subassemblies identified from Gl{A) are subject to further 
evaluation for the selection of preferred direct subassemblies. 

The evaluation of a subassembly is based on the following criteria: 1) 
Stability, 2) Directionality, and 3) Manipulability, which leads to an esti- 
mate of assembly cost involved in local assembly through the identifica- 
tion of the required holding devices and reorientations during assembly. 
Note that the optimal selection of a direct assembly at each stage of 
15 backward assembly planning based on local evaluation may not yield a 
globally optimal plan. As will be shown later, this problem is handled by 
the AO* algorithm with its cost and heuristic functions defined in terms 
of the above criteria. 

5.1 Stability 

20 

The stability of a subassembly can be represented in terms of the internal 
freedom of motion associated with the subassembly. To present a method 
formalizing the internal freedom of motion of a subassembly, let us first 
define the following: A cluster of parts in Sf|A, P*|(Sf|.4). is said to be 
2^ floating if it is connected to the rest of Sf |A only by floating liaisons. 
Pk\(Sf \A) corresponds to a subgraph of Gi(Sf\A) that can be separated 
from G L (Sf\A) by a cut-set consisting only of floating liaisons (denoted 
here as a floating cut-set). A cluster of parts in Sf |.4, P*.|(Sf|A), is 
said to be disconnected, if it has no liaison connected to the rest of 
Sf\A. Pk\{Sf\A) corresponds to a disjoint subgraph in C?/,(Sf|-4). The 
stability of Sf\A can be defined based on a set of floating c lusters of parts, 
P k \(S?\A), and disconnected clusters of parts, P*.|(Sf |-4). included in 
Sf\A. Let us assume that a floating cut-set, f* ? of G'^Sf |.4) decomposes 
Sf\A into Pk\(Sf |A) and Sf|A-P|fc|(Sf|A). The local freedom of motion of 
rjt, LFMifk ), can be defined as the local freedom of motion of i\.|(Sf|A) 
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against Sf\A - Pk\(Sf\A): 

LFM(f t ) = LFM(f%\(Sf\A),St\A- P t \{Sf\A)) 

= - h\(Sf\A)), 

5 

for all /j, lj € r*. 

Note that Pk\(Sf\A) should be chosen in such a way that Sf\A — 
Pk\(Sf\A) includes one of the A-nodes of the assembly .4 that will be 
used for holding Sf\A during the mating of 5f|A with .4 — Sf|-4. Now 
^ let us define the following: 

Definition: Internal Freedom of Motion of Sf\A 

The internal freedom of motion of 5f|A, IFM(Sf |*4), is defined as a 
collection of assembly directions to which Sf\A can be broken apart. 
IFM(Sf\A) can be calculated by the following rules: 

15 1) ? aKSjV) and P t \(Sf\A) =* IFM(Sf\A) = 0; 

2) ?a|(5f|.4) but 3a|(5f|A) =>■ IFM(Sf\A) = U LFU(h), Vr*; 

3) 3P t |(5f|-4) =* IFM(Sf\A) = {±x,±y,±:,± v. ±8.± 4>}- 

As an example, let us consider a simple 2-D assembly shown in Fig. 

20 ? 

Since / 2 and / 3 are floating liaisons, we have that 

ri = {/i,/ 2 }, ?2 — {h,h}, and f 3 = {/ 2 N 3 }. 

Assuming that Sf|.4 is oriented with reference to the assembly pose 
25 of A as shown in Fig. 7, and that Pi, an A-node of .4. is selected for 
grasping of Sf\A, we have 

LFMih) = LPM(P 2 UP 3 ,Pi) = {+x,+c}. 

LFM (r 2 ) = XPM(P 2 ,PiUP 3 ) = {+~}, 

30 LFM{h) = IFM(P 3 ,PiU? 2 ) = {+n+:}. 

Therefore, 

IFM(Sf\A) = {+x,+z}. 

Based on the definition of IFM(S?\A), we can establish the stability 
condition for Sf\A , as follows: 
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1) 5f|A is said to be self-stable, if IFM(Sf\A) is null, or IFM(Sf\A) has 

at most a single translational freedom of motion with a rotational 
freedom of motion only about the axis of translation (i.e. P^\{Sf\A) 
and Sf \A — Pk\(Sf\A) form a single peg-and-hole type mating rela- 
tionship). 

2) is said to be stable with the assistance of holding devices, if 
each P*|(S?|A) or P k \(Sf\A) contains an A- node of the assembly .4. 
This implies that the mating operation of Sf |A can be stabilized and 

10 completed with the assistance of external devices holding P*|(S,- 1-4) 

and Pjfc|(5f|j4) of more than a peg-and-hole type of motion freedom 
against Sf\A - P*|($/|.A). 

3) Otherwise, Sf|A is said to be unstable. 

A stable Sf\A, whether it requires a holding device or not, has one or 
more stable assembly poses, where an assembly pose is represented by 
assembly coordinate aligned with the direction of part stacking against 
gravity. For instance, 5f|A with null IFM h (denoting IFAI after the 
incorporation of necessary holding devices) can have assembly pose of 
Ax, ±y and ±z. 

20 Sf\A with IFM h of {+x,±ip} can have an assembly pose of Ax, re- 

quiring a reorientation of Sf\A to align Ax with the stacking direction 
( against gravity ) . 

Let us now consider the stability associated with an assembly opera- 
tion: 

25 

Definition: Stable Assembly Operation 

The assembly operation between Sf\A and A — Sf\A is said to be sta- 
ble, if Sf\A and A - Sf\A have at least one common stable assembly 
pose defining a stable assembly direction. 

30 Now, the evaluation of Sf | A in terms of stability is based on the sta- 

bility of Sf\A and A - Sf\A , as well as the stability of the assembly 
operation between Sf\A and A — Sf|A: 

1) If either Sf\A is unstable or A - Sf\A is unstable, Sf\A can not be 
selected for a direct subassembly of A. 
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2) When Sf \A and A-Sf\A have no common stable assembly pose, Sf\A 

can not be selected for a direct subassembly of A. 

3) Otherwise we evaluate the assembly cost incurred by the need to 

5 stabilize Sf \A and A—Sf\A as well as the assembly operation between 

Sf\A and A - Sf\A. 

The assembly cost is directly related to the number of holding devices 
required for stabilizing 5f|A and A — SflA, and the necessity of reori- 
enting Sf\A and A — Sf\A for a stable assembly operation. The latter 
will be analyzed in more detail in the next section in relation to the 
directionality in assembly and the determination of best assembly poses. 

5.2 Directionality and Assembly Pose 

The directionality in assembly is another important factor affecting as- 
15 sembly cost. Locally, a stacking operation is considered more cost- 
effective than a non-stacking operation. Globally, a single direction of 
assembly is preferred to multiple directions of assembly. Therefore, the 
evaluation of directionality in assembly should be based on both the lo- 
cal assembly direction between Sf \A and A — 5f |.4, and the uniformity 
20 of assembly directions embedded in Sf |^4 and A — Sf |A 

It should be noted that whether or not the local assembly direction 
between Sf\A and A - Sf\A can be a stacking direction depends on the 
choice of the mating pose (as one of the stable assembly poses common 
to Sf \A and A — Sf |,4). However, the selection of a mating pose between 
25 Sf|j4 and A — Sf |.4 based solely on implementing a stacking operation 
may incur the need to reorient the assembly of Sf\A and .4 — Sf \A. The 
assembly pose of Sf \A or A—Sf\A should be chosen from the set of stable 
assembly poses of Sf\A or A - Sf\A, which may differ from the selected 
mating pose between Sf\A and A — Sf|A. This also incurs the need to 
reorient Sf\A or A — Sf\A, so that Sf\A or A — Sf|-4 is brought into its 
pose. This implies that the determination of an assembly pose and an 
assembly direction should consider the trade-off between maximizing the 
directionality in assembly and minimizing the reorientation of assembly 
pose. In principle, it is desirable to avoid a costly reorientation, unless 
the reorientation is required to allow many local stacking operations in 



22 


the subsequent backward assembly planning, thus justifying the cost of 
reorientation. 

< Cost involved in the Local Mating Operation between Sf\A 
5 and A - Sf\A > 

Let us first introduce the following notational conventions: 

1) tiih = an assembly pose of Sf|A and A — S/|-4 represented with 
reference to the previously determined assembly pose of A, t . 

10 2) {tf}, {i 2 } = a se ^ °f stable assembly poses for Sf\A and .4 — Sf \A. 

3) {tf 2 } = a set of stable assembly poses common to Sf\A and .4 - Sf \A, 
i.e., {t\ 2 } = «} n {<}}. 

We can associate each pair (£f,£ 2 ), where tf £ {tf} and € {t 2 }* with 
15 the relative assembly cost, L, involved in a local mating operation. L can 
be determined based on the required number of reorientations and the 
directionality of mating operations (whether it is a stacking operation or 
a non-stacking operation), as well as the difficulty of handling the related 
subassemblies. 

20 The reorientation of the assembly poses of Sf\A and .4 - S l '\A during 

assembly becomes necessary due to: 

1) The need to transform t\ and/or t s 2 into a mating pose, t * 2 , t j 2 € 12}- 

in the case where t\ ^ t \ 2 or t 2 ^ t\ 2 . 

2 ) The need to transform the selected mating pose, tf 2 , into the assembly 

nr ' 

pose of .4, t*, in the case where t \ 2 ^ t*. 

Table 1 shows the reorientations required for the mating between Sf | -4 
and (A - Sf\A), under various conditions on and t? 2 . The direction- 
ality of the mating operation can be tested by transforming the mating 
3Q directions (between Sf\A and A - Sf\A ) in terms of t* into the mating 
directions in terms of tf 2> where the mating directions in terms of t* are 
identified during the verification of the path existence. 

The relative assembly cost, L, involved in a local mating operation 
can now be calculated for individual (tf,i 2 )i by 


L = [ 'nREO{t \ ) + 7 2 REO(t s 2 ) + lz REO{t\ 2 )]a 0 + - 4 /3 0 , 
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TABLE 1: The Reorientations required for the Local Mating Operation 
between Sf\A and A — Sf\A 


Conditions 

The Required Reorientations 

r t {t; 2 } 

<1 = <2 

t\ € {tf 2 } 

i(<i2 - n 

t\ i {th} 

3(tf ^21^2 — 1 " ^12’ ^12 ~ 4 ^*) 


t\ £ {tf 2 } or 
*2 € {tf 2 } 

2(tf -*• <12 or tf 2 t s v2 J s v2 -* t*) 

t\ £ {«f 2 } and 

<2 t {<12} 


** ^ W2} 

n = t 2 

t\ = t* 

0 

7^ t*, <| £ {^2} 


if £{<;,} 

2(t\ -► t\t% -» /* ) 

*! 

t\ = t* or t s 2 = t* 

1 {t s 2 -*■ t* or t\ t* ) 

7^ t* and t 2 7^ t* 

2(/f -4 t.\t* 2 -> t* ) 


15 

where REO(t\) , REO(t s 2 ), and REO(t\ 2 ) are binary functions of either 1 
(when the reorientation of the corresponding assembly pose is required) 
or 0; ao and /3 q , represent respectively the normal relative assembly cost 
due to a reorientation and a mating motion, 71, 72, 73 and 74 represent the 
20 effect of part manipulability on the relative assembly cost (refer to the 
next section for more detail). The best assembly poses of Sf|.4, t*(Sf |A), 
and A — Sf\A,t*(A — Sf\A) can be determined based on achieving the 
minimum relative assembly cost, L, due to the local mating operation. 

25 < Cost involved in the Assembly of Sf\A and .4 — Sf\A > 

Let us now consider the uniformity of assembly direction embedded 
in Sf\A and A - Sf|A, as a means of globally estimating the relative 
assembly cost, R, involved in the assembly of 5f|A and .4 — 5f|-4. Since 
the exact evaluation of R can only be obtained after a complete assembly 
30 plan is formulated, we indirectly estimate R based on the following two 
major factors contributing to R: 

1) The estimated relative assembly costs, R 0 (Sf\A) and R 0 {A — Sf|-4), 
due to the number of reorientations involved in the assembly of Sf\A 
and A — 5f|A. 
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2) The estimated relative assembly cost, R s (Sf\A) and i? lS (.4 — Sf|.4), 
due to the number of stacking and non-stacking operations involved 
in the assembly of Sf |A and A — Sf\A . 

5 The estimation of R 0 (Sf |A) is based on the directional uniformity of Sf\A 
under the special consideration of the necessary reorientation( flip-over) 
of the bases of Sf\A, as follows: 

Sf\A is said to have m degrees of directionality in £ = x,y , or 2), 

denoted as D^(Sf\A) = m, in the case where there exist m liaisons, 
10 in Gi(Sf\A ) such that +£ or — £ E {LFM(h)} for i = T--m. 

Note that in defining the directionality of Sf \A, x , y, and c are referenced 
in terms of the assembly pose of A , and LFM(lj) can be computed 
either by LFM(l t ;P l ,P 2 ) or LFM{h]P 2 ,Pi ), since LFM(l r , P U P 2 ) = 
— LFM (If, P 2 , Pi). 

15 5f |-4 is said to have the directional uniformity of r in denoted as 

U ( (Sf\A) = r, in the case where Z)^(5f|i4)/Card {/,, /, E G'/.(5f|-4)} = r, 
where r < 1. Sf |-4 has the maximum directional uniformity in £. if Sf \ A 
has the maximum directionality in £. 

Now, let us select U^{Sf\A) from {U^(Sf |A), ^ = x, y, or c} in a de- 
20 creasing order until the accumulation of the selected U^{Sf |-4)"s becomes 
greater than or equal to unity. Define {£*(5/1^4)} as a set containing £’s 
which are associated with the selected U^(Sf\Ays. Then, {^(Sf|.4)} rep- 
resents the minimum number of distinctive assembly directions involved 
in the assembly of Sf|-4, such that Card{f*(Sf|-4)} — 1 provides an es- 
timate of the lower bound of the number of reorientations required for 
the assembly of Sf |-4. A more accurate estimate of the required number 
of reorientations can be obtained by identifying the base flip-overs that 
become necessary when parts are connected to the both side of a base 
such that {LFM(li; Parti, Base), VParti linked to base} includes both 
+£ and — £. To be more specific on this point, let us define the following: 
A base node, rig, of Sf\A , ng(5f|A), is a node of G£(5f|-4) which has 
the degree far greater than the average degree of a node of G'£(-Sf|.4), 
i.e., the degree of n B {S(\A) > k • Average Degree of G L {S-\A). where k 
is a constant such that k » 1. A base node n 5 (5f|,4) is said to have 
the directionality, D nB (Sf\A), where D nB {Sf\A) is the set of independent 
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directions involved in {LFM(li\ P,, Base), V/,-: /,• is linked to tib}, where 
D nB {Sf\A) C {+:r, -x,+y, -y,+z, -z}. In the case where D„b(S?\A) 
contains +£ and — £ (£ = x,y, or z), it is necessary to flip-over the base 


to accomodate +£ and — £. 

Considering the above, R 0 (Sf\A) can be estimated by the following 
equation: 


where 


a i 


h 




Ro(Sf\A) = a 0 [(£ • 6* - 1) + <*], 


2, if +£ and -£ G £>„ B (Sf|A), n B G («b}: 

< A set of base nodes in Gi(Sf\A) 

1 , otherwise 

j 1 , if te{V(S}\A)} 

I 0, otherwise 

f 1, if the selected best pose of Sf|A ^ {^*(5f|A)} 
1 0, otherwise. 


The estimate of relative assembly cost for Sf|A, R(Sf\A), is the addi- 
tion of the estimate of the relative assembly cost due to required reori- 
entations, i? 0 (Sf|A), and the estimate of the relative assembly cost due 
to required part matings, R s (Sf\A ): 

R(Sf\A) = R 0 (Sf\A) + R s (S}\A). 


Since the reorientations considered in the above computation of i? 0 (S/|.4) 
allow part matings to be done by stacking operations, i?*(Sf|A) can be 
estimated simply as 


R $ (Sf\A) = fti-JV, 

where N represents the number of parts included in Sf|.4. 

An alternative way of estimating R(Sf |A) is that R 0 {Sf |-4) is estimated 
based only on the required number of base flip-overs, so as to avoid costly 
reorientations as much as possible. Instead, the assembly of Sf\A allows 
non-stacking operations the cost of which should be incorporated into 
R s (Sf\A), The computation of R 0 (Sf\A) due only to the base flip-overs 
can be computed simply by applying different a ^ to the above equation: 
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1, if +£ and € D nB (Sf\A), n B <E {n B }: 

A set of base nodes in GT(Sf|A) 

0, otherwise 

Whereas, R s (Sf\A ) can be computed by 

R s (Sf\A) = N ■ \f3 m U ( „JSf\A) + Ml - U (m JSf\A))}, 

where 

U Uax = max{ L/f (5f | A), £ = x, y and z}. 

Then, R(Sf\A) is chosen as the lower estimate from the two alterna- 
tives described above. 

Finally, the global estimation of the relative assembly cost, R. due to 
directionality, can be obtained by R = R(Sf\A) + R(A — Sf|A). 

5.3 Manipulability 

A subassembly subject to either a reorientation and/or a translation 
for mating should be easily manipulable by tools or hands. The term 
manipulability of Sf\A is used to quantify the efficiency in orienting Sf\A 
and in handling of Sf\A. The manipulability of Sf\A is closely linked 
to the size, shape and weight of Sf\ A. More specifically, the orientation 
efficiency can be measured based on the symmetry and marked polarity 
in the geometry and weight of SflA, whereas the handling efficiency can 
be measured based on the regularity in the size, weight and shape of 
Sf\A, and the flexibility and fragility of Sf|A, which determine the need 
for special tooling. 

The manipulabilities of Sf\A and A-Sf|A affect the relative assembly 
cost of the local mating operation between Sf\A and A — S?|A, since they 
directly influence the relative assembly cost for the required reorienta- 
tions as well as the mating motion. 

To take this into consideration, in the previous section, the relative 
assembly cost for a reorientation, o;, as well as the relative cost, for a 
mating motion, (3 were determined by multiplying the manipulability 
coefficient, 7, to their nominal values, ao and (3 0 . 



27 


TABLE 2: The Criteria for Measuring Manipulability of a Part or a 


Subassembly 


Orientation Efficiency 

Relative Assembly Cost 

Part tangles, nests or shingles 

5 

Asymmetric part without marked 

5 

polarities of weight or geometry 


Asymmetric part with marked 
polarities of weight or geometry 

3 

Symmetric part 

1 

Part delivered to the assembly station 
with a known orientation 

1 


Handling Efficiency 

Relative Assembly Cost 

Large off center weight potentially 
causing loss of orientation 

5 

Very large parts 

5 

Very small parts 

5 

Fragile 

3 

Flexible 

3 

Irregular shaped part requiring 
special tooling 

3 

Easily handled part with standard tooling 
(tooling can handle more than 1 part) 

1 


The manipulability coefficient of Sf|A, 7(5'f|A), can be determined as 
the ratio between the sum of the score for each criterion of orientation 
and handling efficiency for Sf\A and the sum of the nominal scores for 
each criterion of orientation and handling efficiency. 

30 

6 Selection of Best Subassemblies 
based on AO* Algorithm 

As indicated earlier herein, the selection of Sf\A based solely on the rela- 
tive assembly costs involved in the local mating operation between Sf\A 
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and A — Sf \A may not produce a globally optimal assembly plan. There- 
fore, we adopt the AO* algorithm with a properly defined evaluation 
function to search for a globally optimal or suboptimal plan. 

The search space to which the AO* algorithm is applied can be rep- 
resented by an AND/OR tree. 

The decomposition of an assembly A in backward assembly planning 
implies the expansion of an AND node (representing an assembly A) 
into its OR children representing the alternative decompositions of A, 
{(Sf |A, A — Sf |A), i = 1 ,•••,/}, and its AND grandchildren {Sf \A and 
A — for i = 1, • • • , /} attached to indivisual OR children, as shown 

in Fig. 8. The AO* algorithm searches for an optimal solution tree by 
expanding those AND nodes of the current potential solution tree that 
are open to expansion, and by evaluating the next alternatives based on 
an evaluation function. 

An AND tree is an AND/OR tree every AND node of which has no 
more than one OR child. A potential solution tree is an AND tree having 
the minimum value for the evaluation function at the current stage of 
search, whereas a solution tree is an AND tree with leaves consisting of 
only single parts. 

The Local Cost, c/(nj), associated with an OR node, 7?J, represents the 
relative assembly cost incurred by the local mating operation between 
5f| A and A - Sf\A. c/(n°) can be computed by the weighted sum of the 
following three components: 

1) The relative assembly cost due to the stabilization of 5/|.4 and A — 

5/|-4 by using holding devices and/or reorientations, as described in 
Table 1. 

2) The relative assembly cost due to the reorientations and translations 

required for mating between Sf\A and A — 5/ 1 -4 , as described in 
Table 2 and Table 3. Note that this cost is linked to the directional- 
ities and best assembly poses for 5/1-4 and -4 — 5/|.4. as well as the 
manipulabilities of 5/|A and A — 5/|A. 

3) The relative cost of the special processes assigned to -4, the parent 

node of nf. 
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The Accumulated Cost, c a (T“), associated with an AND tree , T“, rep- 
resents the weighted sum of the following two components: 

1) The sum of c/(n- ) ) for all nj, nj G T“. 

2) The depth of T* defined by the maximum depth of n for all n- 1 , 

n® G T“, where the depth of an OR node is measured in terms of the 
depth among OR nodes without considering AND nodes. 

The Local Heuristic Estimate, h e (n®), associated with an OR node, n®, 
represents an estimate of the optimal relative assembly cost to assemble 
Sf\ A and A — Sf|A, and can be computed by the weighted sum of the 
following components: 

1) The relative assembly cost, R , associated with the directional unifor- 
mity of Sf\A and A — Sf\A, as defined in the previous section. 


2) The relative assembly cost, 5, associated with the internal stability 
of Sf|,4 and .4 - Sj\A : 

S = + x(-4 - Sfl.4)] 

where \(5f|A) and \{A — S^|A) represent the internal stability of 
Sf | A and A — Sf|A, respectively, and are defined by 

the number of floating liaisons in Gl{S‘'\A) 


X{Sf\A) = 
X(A-Sf\A) = 


the average degree of a node in G /. ( Sf \ .4 ) 
the number of floating liaisons in Gl{A — Sf |-4) 


the average degree of a node in Gl(A — Sf |A) 
and S represents the relative assembly cost due to a holding device. 


Then, the evaluation function, e/(T associated with an AND tree 
T* simply becomes 


e f (T t a ) = ^(T?) + r]h e (T, a ), 

where rj adjusts the contribution of h t {T^) to ef(T-‘) in relation with 
Co(T?). 


7 Assembly Planning with DFA Analysis 

As shown in the previous section, the evaluation of the local cost. c/(n®). 
at an OR node, nj, is based on the detailed analysis of c/(n ■’) in terms of 
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the stability, the directionality, the assembly pose and the manipulability 
associated with the assembly of the children of n®, as well as the cost of 
special processes assigned to the parent of n®. 

The result of this analysis at each OR node of the search tree can 
directly be used for the identification of the assemblability of a product 
and for the evaluation of DFA criteria, which can be fed back to the 
designer for proper design evaluation and modification. The assembly 
planner developed here has both the capability of selecting an optimal 
assembly partial order as well as the capability of conducting DFA anal- 
ysis, serving as a powerful tool for automating the DFA evaluation and 
modification cycle in concurrent engineering. 

DFA analysis performed during the process of computing the local 
cost, c/(n®), associated with an OR node, n®, n® ~ {5/|.4..4 — Sf|A}, is 
summarized into the DFA analysis table for n®, as illustrated in Table 3. 

Now, the analysis of DFA for a given product can be accomplished 
based on the DFA tables associated with all of the OR nodes of the so- 
lution tree. 
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Example: The Flashlight Assembly 

The AO* algorithm with the cost and heuristic functions defined in the 
previous section is applied to the flashlight assembly for finding an op- 
timal solution tree and performing DFA analysis. Fig. 9 illustrates first 
several nodes of the AND/OR search tree formed by the AO* algorithm, 
where DFA analysis tables are attached to individual OR nodes. 

At root node(Node 0), which corresponds to the whole product of 
flashlight, there is only one direct subassembly (Node 1). SI or. S'2 im- 
plying that there is only one way to disassemble the product in the first 
step. Node 1 is expanded to it’s AND children, (2,3). Further expansion 
of Node 2 can be done easily because the grouping of B1 and B2 pro- 
vides only one alternative in selecting direct subassembly. At Node 3, 
the abstract liaison graph is shown in Fig. 10. Note that 53 is a group 
node resulted from the special process forest. 

The algorithm to identify direct subassemblies is applied to S2 and 
the result is shown in Table 4. 

Therefore, at Node 3, the system identifies two alternative direct sub- 
assemblies, C and S3. These two alternatives are represented as the OR 
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TABLE 3: A DFA Analysis Table for an OR Node, n°, ~ {5/jA, .4 — 5/|A} 


DFA Analysis Category 

DFA Criteria Details 

Stability 

Total 

relative cost 
due to the 
need to sta- 
bilize 5/|A 
and/or A — 
SflA 

The number of 
holding devices re- 
quired for the sta- 
bilization of 5/| A 
and/or A — 5/|A 

The number of re- 
orientations re- 

quired for the sta- 
bilization of 5/| A 
and/or A — Sf\ A 

The relative 
cost due to t he 
required hold- 
ing devices 

The relative 
cost due to the 
required reori- 
entaion 

Manipulabili 
and Direc- 
tionality 

tyTotal 
relative cost 
involved in 
mating be- 
tween 5f| A 
and A — 
5/| A due to 
manipula- 
bility and 

The manipulabil- 
ity factors 
(Refer to Table 2 
for more details) 

and 

1(A - $f\A) 


The best assem- 
bly poses for 5/| A 
and/or A — 5/|A 



directionality 

The number of re- 
orientations 
required for mat- 
ing between 5/|A 
and/or A — 5/|A 

The relative 
cost due to the 
required reori- 
entations 



The translati- 

nal motion during 
mat- 
ing between Sf\A 
and/or A — Sf\ A 

The relative 
cost 

due to the re- 
quired trans- 
latinal motion 
for mating 

Process 

Total rela- 
tive cost for 
the special 
processes 
assigned to 
A 

The list of special 
processes assigned 
to A 

The relative 
cost of indi- 
vidual special 
processes 


Total cost 
at n° 
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TABLE 4: Valid Cut-Sets or Direct Subassemblies generated from 


10 


Cut-set 

P\A 

A-P\A 

Result 

{h, ho} 

C 

S3+G 

Valid 

S3+G 

C 

Valid 

{^8, ho} 

G 

S3+C 

Failed in the LM test 

S3+C 

G 

Failed in the LM test 

{is, is} 

S3 

C+G 

Valid 

C+G 

S3 

Valid 
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nodes, 5 and 6, in Fig. 9, while the two OR nodes are expanded to their 
AND children, (9,10) and (11,12), respectively. The system then calcu- 
lates the evaluation function at Node 5, e/(Node 5). and the evaluation 
function at Node 6, e/(Node 6), based on the local costs, cy(Node 5) and 
c/(Node 6), and the local heuristic estimates, /i e (Node 5) and h ( (Node 6), 
as follows: 

At Node 5, the system identifies that 

1) C is self-stable, but S3+G needs one holding device due to /*. 

2) The best assembly pose of S3+G needs one reorientation so 
that BC+BL+R is on top of G. The assembly pose of the 
product is initially given in such a way that C is on the right 
side to BD. 

3) C needs one reorientation so that S3+G can be stacked onto 
C. 

4) S3+G+C need to be reoriented after the mating operation of 
S3+G and C. 

5) The manipulability coefficients for S3+G(7i), C(+ 2 ) and 
S3+G+C(7 3 ) are assigned to 0.8, 0.4 and 0.4, respectively, 
based on Table 3. 

Therefore, 

c a (Node 5) = 15 (one holding device) 

+10 x 0.8 + 10 x 0.4 + 10 x 0.4 (3 reorientations) 

+1 x 0.8 (stacking operation) 

+2 (depth of Node 5) 

= 33.8 
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For the calculation of the local heuristic estimate at Node 5, the system 
identifies the following: 

1) The estimates of the relative assembly cost due to the direc- 

^ tional uniformity, R, and the internal stability, S, of C are 

zero, since C is a single part. 

2) S3+C+G consists of 4 parts with the maximum uniform di- 
rectionality of 1 in +x. Therefore, i2o(S3+C+G) = 10 x 1 = 

10 and i? s (S3+C+G) = 1x4 = 4. Asa result, the estimate of 
the relative assembly cost due to the directional uniformity, 

10 R, of S3+C+G is 14. 

3) S3+C+G has 3 floating liaisons and has the average degree 
of node of 2. Therefore, the estimate of the relative assembly 
cost due to the internal stability, 5, of S3+C+G is 22.5 (5 = 

15 x | = 22.5 ), where 15 is used for the relative assembly 

15 cost for a holding device. 

4) The effect of Node 5 on assembly parallelism. can be 
estimated as 6 (w p = 10x3/5 = 6 with 10 assigned as a 
weight). 

As a result, we have that h e (Node 5) = 42.5. Therefore, we have that 

20 

e/(Node 5) = c a (Node 5) + rjh e (NodQ 5) 

= 76.3 (with r) = 1.0). 

At Node 6, the system follows the same steps that are used for calcu- 
lating e/(Node 5), and results are the following: 

25 c a { Node C) = 22.6 

Therefore, we have that 

ey(Node 6) = c a (Node 6) + 7//i e (Node G) 

= 54.6 (with r] = 1.0). 

30 

Finally, comparing e/(Node 5) and e/(Node 6), the system selects 
Node 6 for further expansion. The result of such an expansion is shown 
in Fig. 9. This process is continued until all nodes become single parts. 

Fig. 11 illustrates a process implementing the foregoing, in which 
the merging process is performed (block 110 of Fig. 11). followed by the 
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grouping process (block 120), both of which reduce the number of poten- 
tial direct subassemblies which may be identified. All direct subassem- 
blies axe identified (block 130). Then, an AND-OR tree is constructed 
from a knowledge of all direct subassemblies (block 140). The local cost 
5 of each OR node is computed (block 150), and the tree is "pruned" to 
eliminate the higher cost OR nodes (block 160). Finally, an assembly 
sequence is carried out in accordance with the pruned AND-OR tree. 

Fig. 12 illustrates a system for carrying out the invention, including 
a data base representing the assembly design definitions 180 (such as 
q the liaison graph or map, the assembly constraints and special process 
constraints), a computer 190 which performs the process of Fig. 11 to 
generate a sequence of assembly instructions 200. 

8 Conclusion 
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This specification contributes to automatic assembly planning closer to 
reality by 

1) Developing a backward assembly planner which handles the case where 

an assembly sequence is not the same as the reverse of a disassembly 
sequence. 

2) Improving planning efficiency with the reduction of search space by 

merging and grouping parts based on interconnection feasibility and 
special process precedence constraints. 

3) Establishing assembly process planning by incorporating such special 

processes providing as testing, cleaning, etc, in assembly planning. 

4) Providing subassembly evaluation criteria with a direct connection to 

assembly cost: especially the method for the stability and direction- 
ality of an assembly is presented through which the required number 
of holding devices and reorientations during assembly is identified. 

While the invention has been described in detail by specific reference 
to preferred embodiments thereof, it is understood that variations and 
modifications may be made without departing from the true spirit and 
scope of the invention. 
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BACKWARD ASSEMBLY PLANNING WITH DFA 

ANAYLSIS 

ABSRACT 

5 An assembly planning system that operates based on a recursive de- 

composition of assembly into subassemblies, and analyzes assembly cost 
in terms of stability, directionality, and manipulability to guide the gen- 
eration of preferred assembly plans. The planning in this system in- 
corporates the special processes, such as cleaning, testing, labeling, etc. 
IQ that must occur during the assembly, and handles nonreversible as well 
as reversible assembly tasks through backward assembly planning. In or- 
der to increase the planning efficiency, the system avoids the analysis of 
decompositions that do not correspond to feasible assembly tasks. This 
is achieved by grouping and merging those parts that can not be decom- 
posable at the current stage of backward assembly planning due to the 
15 requirement of special processes and the constraint of interconnection 
feasibility. The invention includes methods of evaluating assembly cost 
in terms of the number of fixtures (or holding devices) 
required for assembly, through the analysis of stability, 
manipulability. All these factors are used in defining 
functions for an AO* search for an optimal plan. 
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NASA Case No. NPO 




(DEFSCHEMA L2 
((INSTANCE LIAISON) 

(CLASS FLOATING) 

(MATING PARTS P2 P3) 
(INTERCONNECTION ATTACH) 
(FEATURES(F1 .P2 F1.P3(+X -X))) 
(CONNECTION TYPE PLACE ON))) 
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(DEFSCHEMA P3 (DEFSCHEMA F2.P14 

((INSTANCE PART) ((INSTANCE FEATURE) 

(CLASS GENERAL) (FEATURE OF P3) 

(FEATURE-LIST FI .PI 4 F2.P14) (TYPE CONNECTION) 
(COORD(7 0 10 0 90 0)) (SHAPE(CYLINDER 6 1.8)) 

(VOLU M E((CYLI N DER 14 1) (LOCATION^ 0 0 0 0 0)))) 

0 0 0 0 0 0 ))) 


FIG. 2 
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